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We investigate the influence of static correlations beyond the pair level on the dynamics of selected 
model glass-formers. We compare the pair structure, angular distribution functions, and statistics 
of Voronoi polyhedra of two well-known Lennard- Jones mixtures as well as of the corresponding 
Weeks-Chandler- Andersen variants, in which the attractive part of the potential is truncated. By 
means of the Voronoi construction we identify the atomic arrangements corresponding to the locally 
preferred structures of the models. We find that the growth of domains formed by interconnected 
locally preferred structures signals the onset of the slow dynamics regime and allows to rationalize 
the different dynamic behaviors of the models. At low temperature, the spatial extension of the 
structurally correlated domains, evaluated at fixed relaxation time, increases with the fragility of 
the models and is systematically reduced by truncating the attractions. In view of these results, 
proper inclusion of many-body static correlations in theories of the glass transition appears crucial 
for the description of the dynamics of fragile glass-formers. 

PACS numbers: 61.43.Fs, 61.20.Lc, 64.70.Pf, 61.20.Ja 



I. INTRODUCTION 

The question of whether the structure of a glass differs 
from that of the corresponding liquid is often rhetorically 
posed vifithin the glass community. In fact, very small dif- 
ferences are observed in the static structure factor of a 
viscous liquid as it approaches the glass transition tem- 
perature Tg. By contrast, the viscosity and structural 
relaxation times increase by several orders of magnitude 
upon supercooling and the motion of the molecules in 
the liquid becomes increasingly cooperative and spatially 
heterogeneous. 

At first glance, the small structural changes discernible 
at the level of pair correlations appear insufficient to 
explain the dramatic slowing down of the liquid and 
the non-trivial spatial correlations of the dynamics. A 
preliminary indication that this may not necessarily be 
the case already comes from the Mode Coupling theory 
(MCT) of the glass transition [T]. The MCT predictions 
for the dynamic correlation functions are based uniquely 
on structural information, almost invariably the static 
structure factors of the liquid. Numerical solutions of 
the MCT equations for model liquids show that small 
variations of the pair correlation functions, which de- 
velop upon lowering the temperature, can produce signif- 
icant effects on the dynamics and eventually lead to an 
ergodic-non-ergodic transition at some critical tempera- 
ture Tmct p] • The generic predictions of the theory ac- 
count rather well for experimental and numerical findings 
in weakly supercooled liquids [HIS] (i.e., for T > Tmct)- 
The break-down of MCT at lower temperatures {Tg < 
T < Tmct), where the actual system remains effectively 
ergodic, most likely reflects the mean-field character of 
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the theory |3HS] and its inability to describe activated 
transitions between metastable glassy states [7] . Another 
delicate aspect that may affect the outcome of MCT is 
the exclusion of many-body correlations |5]. Three-body 
static correlations have been shown [9] to impact the 
MCT solutions for a model of silica [10], but not the 
ones for the prototypical Kob- Andersen model [TT]. In- 
terestingly, a study of a schematic version of the gener- 
alized MCT [T2Ul4] . which allows proper description of 
many-body dynamic correlations, shows that the ideal 
transition at Tmct can be delayed by retaining higher 
order density correlations in the MCT equations |T3]. 

The importance of high order static correlations, hid- 
den in the amorphous structure of the liquid, is particu- 
larly emphasized by frustration-based approaches to the 
glass transition [TSHTS] . According to these theories, the 
phenomenology of glass formation arises from the com- 
petition between the growth of slow, correlated domains, 
characterized by some preferred local order, and frus- 
tration, which prevents these domains from percolating 
through the liquid. Despite some disagreement on the 
interpretation of the role of frustration [19l [20], these 
models indicate medium range order and structural cor- 
relations beyond the pair level as key features for under- 
standing the dynamic behavior of glass-forming systems. 

Computer simulations of several model glassy sys- 
tems pTI - l?f] and experiments on dense colloidal suspen- 
sions [281 129] provide evidence for the existence of do- 
mains formed by preferred local structures and for their 
influence on the dynamics. Similar observations, albeit 
without explicit reference to the dynamics, emerge from 
recent ah initio simulations and experiments on metallic 
glasses [3DJ [31]. Furthermore, high order static corre- 
lations, named "point-to-set" correlations [32 , have re- 
cently been revealed by simulations under amorphous 
boundary conditions and have been found to grow by de- 
creasing temperature in a model supercooled liquid [33j . 
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FIG. 1. Main panel: pair potentials between particles of 
species 1 used in this study: LJ (dashed line), WCAS (solid 
line), and WCA (dotted line). Inset: force between particles 
of species 1 for WCAS (solid line) and WCA (dotted line) 
potentials. In both panels, the dotted vertical line marks the 
distance r = 2^^^ corresponding to the minimum of the LJ 
potential. 



In spite of these advances, there is still no general consen- 
sus on the connection between the structure and dynam- 
ics in supercooled liquids. In particular, dynamic facili- 
tation models |34j provide an alternative and physically 
appealing description of the glassy dynamics in terms of 
purely kinetic constraints. 



A clear-cut procedure to test the influence of many- 
body static correlations on the dynamics of glass-forming 
liquids emerges from recent work of Berthier and Tar- 
jus [Hiss]. These authors compared the pair structure 
and dynamics of two model glassy systems: the Kob- 
Andersen (KA) binary Lennard- Jones (LJ) mixture [11] 
and its Weeks-Chandler- Andersen (WCA) variant 
in which the attractive part of the pair potential is trun- 
cated. Berthier and Tarjus found that, at fixed tempera- 
ture and for sufficiently large density, the pair structure 
of the two models is almost identical, while the struc- 
tural relaxation times can differ by orders of magnitude. 
Thus, direct comparison of LJ and WCA models offers 
an ideal benchmark to test the existence and the influ- 
ence of static correlations beyond the pair level. Build- 
ing on prior knowledge on the preferred local order of 
L J mixtures |26| , we will consider here structural indica- 
tors of increasing complexity — ranging from pair correla- 
tions, through angular distribution functions, to Voronoi 
tessellation — for selected LJ and WCA liquids and per- 
form a crucial numerical experiment on the influence of 
structure on the dynamics of the models. 



II. METHODS 

We will consider two well-known models of glass- 
forming liquids: the Kob-Andersen binary mixture [llj 
and the Wahnstrom binary mixture (WAHN) |37j. In 
the original models [11] 137], named herein KA-LJ and 
WAHN-LJ, particles interact through the LJ potential 



Uapir) = 4e„/3 



(1) 



where a, f5 = 1,2 are species indices. The values of the 
parameters cTq^ and e^^ are ai2 — O.SfJn, 0-22 = O.SScrii, 
£12 = 1.5eii, and £22 = O.Sen for the KA mixture, and 
CT12 = 0.916crii, a-22 = 0.833ctii, and £22 = £12 = en 
for the WAHN mixture. The chemical compositions and 
mass rations are Xi = 1 — X2 = 0.8, mi/m2 — 1 (for 
the KA mixture) and xi = X2 = 0.5, 7711/7712 = 2 (for 
the WAHN mixture). The potentials are cut and shifted 
by a quadratic term |38j at 2.5fTQ^ and 2.5tTii in KA 
and WAHN mixtures, r espectively. In the following, we 
will use (Til, £11, and miali/en as units of distance, 
energy, and time, respectively. 

In addition, we study the corresponding WCA variants 
of the above mixtures. In the WCA models [Ml 131], the 
interaction parameters CTq^ and £q^ and chemical com- 
positions are unchanged but each of the pair potentials 
Uapir) is truncated and shifted so that the value at the 
minimum is zero [31] . The WCA truncation of the attrac- 
tive part of the potential is well-known from liquid state 
theories 39]. We found, however, that this procedure 
leads to poor energy conservation during the long molec- 
ular dynamics simulations in the supercooled regime. To 
circumvent this problem, we employ here a smooth cut 
off scheme with cubic interpolation |40j to ensure conti- 
nuity up to the second derivative of the potentials at the 
minimum — 2^/^(Tq^ of the LJ potential. Explicitly, 
the WCA smoothed potentials (WCAS) read 



r -Uapir) + Aafi r < Gap 

j(r) = <( Bapirc-r)^ a^p < r 




(2) 



where Aap and Bap are determined to ensure continu- 
ity at r = a and r — re- The parameters aap are 
adjusted for each pair a-/? so that — 2^^^aap, and 
read an = 1.0269, 012 = 0.8215, 022 = 0.9038 for the 
KA-WCAS mixture, and an = 1.0269, ai2 = 0.9473, 
022 = 0.8555 for the WAHN- WCAS mixture. A com- 
parison between LJ, WCA, and WCAS potentials for 1-1 
pairs is shown in Fig. |T] In contrast to the WCA poten- 
tial, the derivative of the force of the WCAS potential is 
continuous at r^. In the inset, we highlight the difference 
between the WCA and the WCAS potentials in the r ~ rj, 
region. As it will be clear in the following, this modifica- 
tion introduces some small differences in the thermody- 
namic and dynamic properties, but does not qualitatively 
alter the comparison with LJ models. In the following, we 
will mostly concentrate on the WCAS models and report 
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FIG. 2. Radial distribution functions gi2{r) for (a) the Kob- 
Andersen mixture and (b) the Wahnstrom mixture. The state 
points considered are (a) p — 1.2, T = 0.5 and (b) p — 1.297, 
T = 0.6. In both panels, solid, dashed, and dotted lines 
indicate results for the LJ, WCAS, and WCA models, respec- 
tively. Error bars are smaller than the widths of the lines. 



selected results for the original WCA models. All stud- 
ied systems are composed of 1000 particles in a cubic box 
with periodic boundary conditions. Molecular dynamics 
simulations are performed in the NVT ensemble using 
the Nosc-Poincare thermostat [H] with a mass parame- 
ter Q=5.0. The number density of the KA mixtures is 
p = 1.2, while that of Wahnstrom mixtures is p = 1.297. 
For the LJ and WCAS models, static and dynamic prop- 
erties are averaged over up to six independent thermal 
histories. We find that the KA-WCAS mixture crystal- 
lizes more easily than the other systems [42j . A similar 
tendency to crystallize has been reported in Ref. [13] for 
the KA-WCA model. Only the non-crystallizing samples 
are retained to perform the averages. 



III. RESULTS 

A. Two-body and three-body static correlations 

To start the discussion, we analyze the pair structure 
of the present models. Figure [2] displays the radial dis- 
tribution function §12(1^) between unlike species in the 
KA mixtures (top panel) and WAHN mixtures (bottom 
panel) . For each type of mixture, the results obtained are 
shown at a common temperature, representative of the 
slow dynamics regime of the LJ models. The pair struc- 
ture appears essentially unaffected by the truncation of 
attractions, thus confirming the observations of recent 
simulation works [HI Hi] ■ This result holds for both WCA 
and WCAS models. Only a close inspection of the figures 
reveals that the first minima of gi2{r) are slightly deeper 
in the LJ models, suggesting that the latter systems are 
effectively more supercooled. A similar effect is visible in 



FIG. 3. Angular distribution functions Di2i(f) for (a) 
Kob-Andersen mixtures and (b) Wahnstrom mixtures. State 
points and lines are the same as in Fig. [2] Error bars are 
smaller than the widths of the lines. 



the radial distribution functions reported by Pedersen et 
al. [U for KA-LJ and KA-WCA mixtures. 

At the temperatures considered in Fig. [2j the struc- 
tural relaxation times of the LJ and WCA models differ 
by almost two orders of magnitude, as is evident from 
Ref. '5] and Figs. |4] and [5] (discussed in further detail 
below). Given the small differences observed in the pair 
structure, it is natural to ask whether this large variation 
in the dynamic properties is due to static correlations of 
higher order. In a first attempt to go beyond pair correla- 
tions, we calculate the bond-angle distribution functions 
Da/j-yiP) between triplets of neighboring particles of par- 
ticles of species a, /3, and 7, where (3 is the species of the 
central particle. Figure [3] shows the angular distribution 
functions Di2i{9) for the same state points considered 
in Fig. [2] Angular correlations reveal more clearly the 
structural differences between the LJ and WCA mod- 
els. In KA mixtures, the sharp peaks in Di2i{9) around 

70° and the broad peak in the range 120° - 140° reflect 
local arrangements corresponding to distorted twisted 
bi-capped prisms of large particles (species 1) centered 
around small particles (species 2) [21]. A comparison of 
the LJ and WCA data sets thus reveals that the KA-LJ 
mixture has a more pronounced local ordering than the 
KA-WCA mixture at the selected thermodynamic state. 
As in the case of gi2{'r), the difference between the WCA 
and WCAS models is negligible for this state point. A 
similar effect is visible for the WAHN mixture: the peaks 
in Di2i{9), located around 63°, 116°, and 180°, are sig- 
natures of local icosahedral ordering, which appears more 
pronounced in the original LJ model than in the WCAS 
variant. Similar conclusions can be drawn from an analy- 
sis of the other angular distribution functions (not shown 
here) and are corroborated by an inspection of data at 
even lower temperature. We conclude that the structure 
of LJ and WCA systems differ more evidently at the level 
of three-body static correlations and that the increase 
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FIG. 4. (a) Structural relaxation times r as a function of 1/T 
for Kob- Andersen mixtures: LJ (filled circles), WCAS (open 
circles), and WCA (open triangles) models. The wave- vector 
considered for the calculation of r is fc = 7.0. Fits to the 
modified VFT equation (Eq. ^) are shown as solid lines, (b) 
Average number A^'lps of particles in LPS domains formed 
by (0,2,8)-polyhedra. Symbols have the same meaning as in 
(a), (c) Average fraction of particles of species 2 at the center 
of (0,2,8) polyhedra as a function of l/T. Symbols have the 
same meaning as in (a). 



of local ordering upon switching on attractions correlate 
qualitatively with the increase of relaxation times. 



B. Locally preferred structures 

To render the connection between the local struc- 
ture and dynamics explicit, we analyze the statistics 
of Voronoi polyhedra as a function of temperature. 
The Voronoi tessellation implicitly entails more complex 
static correlations (although, of course, it cannot be ex- 
pressed as a multi-particle correlation function) and re- 
veals the details of the particles' arrangements within the 
first coordination shell. Inspection of the spatial persis- 
tence of a given local structure provides information on 
extended structural correlations, i.e., medium range or- 
der. The protocol adopted here is the same as that used 
in a prior investigation of the local structure of binary 
LJ mixtures at constant pressure |26j . We monitor the 
temperature dependence of the fraction of Voronoi poly- 
hedra with a given signature (713,^4,...), where is 
the number of faces of the polyhedron with a given num- 
ber i of vertices. We identify the locally preferred struc- 
ture (LPS) of the liquid as the geometrical arrangement 
corresponding to the most frequent Voronoi polyhedron 
around particles of species 2 observed in the samples at 



FIG. 5. Same as Fig. |4] but for Wahnstrom mixtures. For 
these systems, the LPS corresponds to (0,0,12) polyhedra. 



low temperature. This choice is based on the observation 
that in binary L J mixtures it is easier to characterize local 
order around small particles H5] . This is also consis- 
tent with a previous study of the KA-LJ mixture, which 
focused on the coordination polyhedra of large particles 
around the small ones 05] . The most frequent signatures 
of Voronoi polyhedra around particles of species 2 and 1 
are reported in Table [l] and [TTl respectively. We find that 
the typical Voronoi polyhedra observed around particles 
of species 1 [such as (0,2,8,4) or (0,1,10,4) polyhedra] do 
not display evident symmetries and lack a clear structural 
identification. Furthermore, the corresponding percent- 
ages do not increase as sharply upon decreasing temper- 
ature as those calculated for Voronoi polyhedra around 
particles of species 2. In the following, we will therefore 
base our analysis on the local structures observed around 
this latter type of particles. Understanding the nature of 
local order around large particles remains an open issue, 
which may need more refined methods to detect short 
and medium range order. 

Our results confirm the observations of Ref . [^B] at con- 
stant pressure: in the low temperature regime, (0,2,8) 
and (0,0,12) polyhedra around particles of species 2 con- 
stitute the dominant signatures in KA-LJ and WAHN-LJ 
mixtures, respectively. Thus, we identify the LPS of KA- 
LJ and WAHN-LJ as twisted bi-capped square prisms 
and icosahedra, respectively |26j . The identification is 
consistent with previous investigations on KA-LJ clus- 
ters [37] and with a recent simulation study on the bulk 
WAHN-LJ mixture f48| . By applying the same procedure 
to the WCA and WCAS models, we find that the locally 
preferred structures remain the same as in the original 
LJ models. We observe, however, a systematic reduction 
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TABLE I. Most frequent signatures of Voronoi polyhedra around particles of species 2 from instantaneous configurations and 
local minima of potential energy surface. Percentages are calculated with respect to the number of particles of species 2. 
The low temperature data set (T — Ti) corresponds to the lowest available temperatures: T = 0.435 (KA-LJ), T = 0.285 
(KA-WCAS), r = 0.560 (WAHN-LJ), and T = 0.343 (WAHN-WCAS). The high temperature data set (T « T*) corresponds 
to temperatures close to the crossover temperature T*: T = 0.983 (KA-LJ), T = 0.627 (KA-WCAS), T = 1.072 (WAHN-LJ), 
and T = 0.598 (WAHN-WCAS). 
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TABLE II. Same as Table |T] but for Voronoi polyhedra around particles of species 1. Percentages are calculated with respect 
to the number of particles of species 1. 
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of the fraction of LPS upon truncating the attractions. 
This effect will be discussed in further detail below. 

As a general rule, the fraction /lps of particles of 
species 2 at the center of a LPS increases with decreasing 
temperature . The growth of /lps reflects the forma- 
tion of slow, long-lived clusters of neighboring LPS [26] . 
A calculation of the self intermediate scattering functions 
filtered according to the pertinent Voronoi polyhedra, 
shows that the typical relaxation times of particles at 
the center of LPSs are up to 10 times larger than those 
outside LPSs . In the following, we will refer to these 
clusters as "LPS domains", which are defined as groups 
of particles sitting either at the center or on the vertices 
of face-sharing polyhedra with the signature of the LPS. 
The average number of particles forming an LPS domain 
will be denoted by A'lps i which is a measure of the spa- 
tial extension over which the liquid adopts the same pre- 
ferred local structure [49^. In Figs. |4|^b) and [sjb) we 
show TVlps as a function of l/T for the KA and WAHN 
models, respectively. To facilitate a comparison with pre- 
vious work [26], we include the temperature dependence 
of /lps in Figs, ^c) and [Sjjc). Note that, while /lps 
is evaluated with respect to particles of species 2, both 
species of particles contribute to the size iVLPs of LPS do- 
mains. Both A^LPS and /lps increase in a similar fashion 
as T decreases, although with slightly different functional 
forms. The growth of LPS domains is particularly dra- 
matic in WAHN mixtures, which develop a strong icosa- 
hedral order upon supercooling. By contrast, the size 
of the domains formed by prismatic structures in KA 
models is relatively small (20-30 particles). Nonetheless, 
the structural evolution in all the systems studied follow 
qualitatively similar patterns. 



C. Connection between structure and dynamics 

To illustrate the connection to the dynamics of the 
models, we now study the temperature dependence of 
the structural relaxation times r. The latter are defined 
by the condition Fs{k = 7, t) — 1/e, where Fs{k,t) is 
the self intermediate scattering function averaged over 
all particles. The relaxation times have been fitted by 
the following modified Vogel-Fulcher-Tammann (VFT) 
equation [25] 



r(T) 



> exp 



, exp 



Eoo/T] 



K{T/To - 1) 



where 



, exp 



Eoo/T* - 



1 



K{T*/Tq - 1) 



(3) 



(4) 



Equation ([s]) ensures a smooth crossover around T* be- 
tween the Arrhenius law at high T and the VFT equation 
at low T, accounting for the super- Arrhenius dependence 
of the relaxation times. Figures Qa) and [sf^a) display r 



as a function of l/T for the KA and WAHN mixtures, re- 
spectively, together with the corresponding fits to Eq. (|3| . 
Figure |4|^a) also includes results for the KA-WCA mix- 
ture obtained using the original cut and shift at the min- 
imum of the potentials, as in previous works [51 [551 [5S] . 
The latter data set is in good agreement with the re- 
sults obtained in Refs. [SJ [3S] for a similar wave-vector 
(k = 7.21). At sufficiently low temperature, however, 
non-negligible deviations appear with respect to the KA- 
WCAS mixture. This discrepancy may be attributed to 
the modification induced by the smooth cut off employed 
in this work. The comparison between the LJ and WCAS 
models, however, remains qualitatively unaffected and 
confirms the general conclusions of Ref. [3 135] . 



A comparison of the LPS analysis and relaxation times 
data reveals a striking relationship between structure and 
dynamics. The increase of r below the crossover temper- 
ature T* correlates to the increase in size of the LPS 
domains. This connection is particularly evident for the 
two WAHN mixtures, in which the increase of icosahe- 
dral order upon decreasing T is very sharp. Our results 
reveal that the large difference in the dynamic behavior 
between the LJ and WCA models reflects different stages 
of the evolution of the local structure on the way to glass 
formation — an effect that is barely visible at the level of 
pair correlations. It should also be noted that the size of 
the LPS domains at low temperature is slightly larger in 
the WCA than in the WCAS models, which is consistent 
with the discrepancy in relaxation times observed above. 



A remark on the nature of local order in these models is 
in order. It has been shown recently that the WAHN-LJ 
mixture can phase-separate and then partially crystal- 
lize in a complex crystal structure that accommodates 
distorted icosahedral geometries [55'. The LPS observed 
in the liquid is thus analogous to the typical local struc- 
ture of the underlying crystal, which is at odds with the 
paradigm of the frustration- limited domains theory [16) . 
The situation is less clear in the case of the KA-LJ mix- 
ture, for which an unambiguous identification of the crys- 
talline phase is lacking. Previous studies [SO] [51] have 
shown that, for chemical compositions close to the one of 
the original model, stable crystals either have CsCl sym- 
metry or are composed of a mixture of fee and hep struc- 
tures of large particles. Interestingly, we found that the 
KA-WCAS model can crystallize into a fee lattice of large 
particles. In this model crystallization is associated with 
a sudden drop in the fraction of (0,2,8) polyhedra and 
a rapid increase of (0,4,4,6) polyhedra centered around 
particles of species 1. This reveals a potential mismatch 
between the LPS of the liquid and the typical local struc- 
ture of the crystal. The question of whether the locally 
preferred structure should coincide with the structure of 
the underlying crystal [HI [20] certainly deserves further 
investigations. 
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FIG. 6. Average fraction /lps of LPS in local minima (open 
circles, left axis) and average potential energy Um of local 
minima (filled circles, right axis) as a function of T for Kob- 
Andersen mixtures: (a) LJ model and (b) WCAS model. 
The dynamic crossover temperatures T* obtained from fits 
to Eq. (|3| are indicated as vertical dotted lines. 
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FIG. 7. Same as Fig. [6] but for Wahnstrom mixtures. 



D. Connections between structure, potential 
energy landscape and fragility 

The connection between the growth of LPS domains 
and slow dynamics is further corroborated by the anal- 
ysis of the potential energy landscape (PEL). It is well- 
known that the appearance of super- Arrhenius behavior 
and non-exponential relaxation around the so-called on- 
set temperature Tq coincides with a sharp change in the 
properties of the local minima of the PEL explored by 
system [S2]. In fact, the average energy u„i{T) of local 
minima remains almost constant at high T and starts 
decreasing rapidly below Tq ■ Figures [6] and [7] display 
Um and the fraction /^pg of particles of species 2 at 
the center of the LPSs evaluated for local minima of 



the PEL as a function of temperature for the KA and 
WAHN mixtures, respectively. We note that the values 
of T* obtained from Eq. [3] are only slightly larger than 
the onset temperatures estimated from the appearance 
of two-step, non-exponential relaxation in the dynamic 
correlation functions, and are consistent, at least for KA 
mixtures, with the values of Tq reported in Ref. |53) . 

Strikingly, the onset of slow dynamics, indicated by 
the drop in Um, always correlates to a sharp increase of 
/lps- We remark that the percentages of other, less fre- 
quent signatures of Voronoi polyhedra do not increase 
as sharply across T* , although some of them do display 
some change upon decreasing T (see Table |l]). We thus 
attribute the onset of slow dynamics in the models stud- 
ied herein to the growth of structural correlations. We 
also note that landscape approaches based on high-order 
stationary points of the potential energy surface |54j may 
provide a complementary view on our results. We found 
in fact that particles at the center of LPSs participate 
less to unstable modes of saddle points [55H57] . Estab- 
lishing a clear relationship between the growth of LPSs 
and the disappearance of unstable directions in the land- 
scape is an interesting open issue that is left for future 
investigation. 

The present findings ostensibly indicate that the drop 
in Um is connected to the formation of energetically fa- 
vored structures. This is consistent with the observation 
that the potential energy associated with the LPS is typ- 
ically lower than that of other structures |55j . In general, 
however, specific local structures may be favored also for 
non-energetic reasons, such as more efficient packing or 
non-trivial entropic effects (e.g., favorable arrangements 
of interlocking LPS) . An interesting example of this com- 
petition is provided by a model of a NiY alloy based 
on LJ interactions |3S^, for which the LPS — a capped 
trigonal prismatic structure — corresponds to a Voronoi 
polyhedron having the smallest volume but not the low- 
est potential energy [55!. Successful attempts to account 
for this complex interplay at the mean-field level in the 
one-component LJ liquid [58j and in a soft-sphere mix- 
ture [Sni [50] must also be acknowledged. 

To set forth the present results in a more compact fash- 
ion, we now plot r as a function of A'^lps (see Fig. |8|, 
thereby making the temperature dependence implicit. 
This representation of the data allows us to illustrate 
more clearly several system-specific aspects of the rela- 
tionship between the structure and dynamics. As ex- 
pected, we find that the increase of r correlates to that 
of A^LPS- However, the spatial extension of LPS domains 
at fixed relaxation times increases systematically with the 
fragility K of the model, which is estimated from fitting 
the relaxation times to Eq. [3] In the case of the more 
fragile WAHN mixtures, the increase of relaxation times 
is evidently dictated by the growth of LPS domains. Over 
the same range of r, the less fragile KA mixtures show 
a weaker structural change, indicating that other effects, 
such as dynamic facilitation, may also be playing an im- 
portant role. The overall trend of variation is consistent 
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Nlps 

FIG. 8. Relaxation times r as a function of the average num- 
ber A^LPS of particles in the LPS domains. The corresponding 
fragility indices K of the models, obtained from fits to the 
modified VFT equation [Eq. ([3|], are also indicated. 

with the correlation between the fragility and thermal 
rate of growth of LPSs put forth in Ref. [35] and suggests 
that the impact of static correlations on the dynamics 
should be more pronounced the more fragile the liquid. 
We also note that the inclusion of attractions tends to in- 
crease A'lps a-t fixed relaxation times. This stabilization 
effect is in qualitative agreement with recent observations 
on LJ and WCA fluids close to the triple point [61] . 



IV. CONCLUSIONS 

We have performed a crucial test on the dynamic role 
of static correlations in glass-forming liquids by compar- 
ing two well-studied LJ mixtures and their corresponding 
purely repulsive variants. Truncation of the attractive 
part of the potential considerably shifts the glass tran- 



sition to lower temperatures and reduces the fragility of 
the liquid, but does not alter the pair structure signif- 
icantly [35]. These phenomena have been explained by 
resorting to indicators revealing more complex structural 
correlations. Building on prior work ,26j , we have identi- 
fied correlated domains formed by locally preferred struc- 
tures. We have found that the growth, by decreasing 
the temperature, of LPS domains is tightly connected to 
the onset of the slow dynamics regime. Furthermore, an 
analysis of LPS domains has allowed us to clearly dis- 
tinguish the different dynamic behaviors of the LJ and 
WCAS models in terms of their structure. In retrospect, 
these results suggest that even small differences that are 
discernible at the level of pair correlations, may reflect 
substantially different stages of the structural evolution 
of a supercooled liquid and be associated with very differ- 
ent dynamic regimes. A solution of the MCT equations 
for the dynamic correlation functions — using two-body 
static correlations as input — only partially accounts for 
the different dynamic behaviors of the L J and WCA mod- 
els [3 112]. Thus, proper inclusion of many-body static 
correlations in theories of the glass transition based on 
structural information seems crucial for a correct descrip- 
tion of the dynamics in fragile glass-formers. This sug- 
gests revisiting and extending previous attempts [9] along 
these lines based on mode-coupling theory. Investigations 
of high order static correlations extracted from simula- 
tions under amorphous boundary conditions [33j . using 
patch repetition [63 or order mining methods [64' , as well 
as implementation of alternative methods for LPS deter- 
mination [551 IS3 1 may provide further clues to improve 
our theoretical understanding of the glass transition. 
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